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Abstract 

We solve the geometry of the closed string field theory five-point vertex. Our solution is calcu- 
lated in terms of quadratic Strebel differentials which are found numerically all over the relevant 
subspace of the moduli space of spheres with five punctures. Part of the boundary of the reduced 
moduli space is described in terms of an algebraic curve, while the remaining part has to be eval- 
uated numerically. We use this data to compute the contact term of five tachyons and estimate 
its uncertainty to be of about 0.1%. To put to a test the theory and the computations done, we 
calculate the contact term of five dilatons. In agreement with the dilaton theorem, it is found to 
cancel the term obtained from the tree level Feynman diagrams built with three- and four- vertices. 
This cancellation, achieved with a precision of about 0.1%, is within the estimated margin error 
on the contact term and is therefore a very good evidence that our computations are reliable. 
The techniques and numerical algorithm developed in this paper make it possible to compute the 
contact amplitude of any five off-shell closed bosonic string states. 
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1 Introduction 

The object of this paper is the explicit computation of the quintic term of the action of closed bosonic 
string field theory (CSFT). This action, formally constructed in [H El [3l H] , is nonpolynomial. This 
is in contrast with Witten's string field theory [5] which is cubic, i.e. Feynman diagrams constructed 
with three- vertices are enough to cover (exactly once) the whole moduli space of Riemann surfaces 
with N punctures on the boundary. In CSFT it is not possible to do so, The Feynman diagrams 
with closed three-vertices do not suffice to construct all spheres with four punctures, one has to 
introduce, in the action, a contact term of order four in the string field to account for the remaining 
four-punctured spheres. But this is not yet enough, one has to introduce also a five-vertex because 
the Feynman diagrams with three- and four- vertices do not cover the whole moduli space of spheres 
with five punctures. And so on, one must put in the action, contact terms of all orders. In this 
paper we will discuss only the classical action, but if we were to consider diagrams with loops, we 
would face a similar problem; namely the vertices of the classical action are not enough to generate 
all Riemann surfaces of genus one; and so on, one must introduce terms at all genii. 

Being able to explicitly compute the whole action would be immensely useful in understanding 
nonperturbative physics of closed strings. The hot problem that we have in mind is, in particular, 
to understand if closed bosonic string theory has a stable vacuum, and what would the theory look 
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like in this vacuum. The approach that we are taking, is by order truncation of the CSFT action, 
i.e. we truncate it to a polynomial of a given order in the string field. The simplest nontrivial 
truncation is done by keeping the quadratic and cubic terms only. In [6j, Kostelecky and Samuel 
did precisely that; they then truncated the string field itself by keeping only the tachyon and all 
the massless fields. In this approximation they found a locally stable vacuum in which the tachyon 
had a positive expectation value. 

The next order of approximation is to keep the quartic term as well. The computation of this 
term is already seriously complicated. The reason lies in the fact that, in order to compute a 
contact amplitude with external states, one has to integrate a certain correlator over a region 
of the moduli space of spheres with N punctures (we call this region the reduced moduli space, 
it corresponds to all spheres that cannot be constructed with Feynman diagrams), which has real 
dimensionality 2(A^ — 3). The expression of the correlator at a point in the reduced moduli space, 
depends on the geometry of the vertex at this point, which in turn is given by the solution of a 
minimal area problem [4]. For the cubic vertex = 3, there is no moduli space to integrate over 
and it is thus easy to calculate. But for the quartic vertex, there is a two-dimensional reduced 
moduli space. This vertex was solved numerically in [7J. The solution given there consisted of 
the boundary of the reduced moduli space in the complex plane, and everywhere in this region 
the geometry of the vertex was expressed with a quadratic differential (see O [9] for details on 
quadratic differentials) given in terms of a complex parameter a(^, ^) depending on the coordinates 
on the reduced moduli space. The solution was explicitly given by a reasonably short fit that can 
be copied from the paper and used to compute amplitudes with an accuracy of about 0.1%. 

The results of [7] were checked by Yang and Zwiebach in [101 lllj . For this, they verified that 
the quartic term in the effective potential of some marginal fields is seen to vanish as one increases 
the truncation level of the string field. A similar analysis was made for the effective potential of 
the dilaton. In particular, they found that the contact quartic term cancels the terms from cubic 
vertices with a precision of about 0.2%. 

In [121 113j . Yang and Zwiebach went on to address the question raised before, whether closed 
bosonic string theory has a stable vacuum. They started by realizing that the tachyon condensate 
must drive the zero-momentum ghost dilaton. Indeed this state with a peculiar ghost structure, 
given by 

= (cic_i - cic„i) |0) (1.1) 

has to be included in the condensate as soon as one considers the quartic vertexj^ This stems 
from the fact that the antighost insertion in the correlator, can make the ghost numbers work so 
that amplitudes with a single dilaton, for example three tachyons and a dilaton, are nonzero. They 
truncated the string field to level four, including the tachyon (level zero), the dilaton (level two) and 
four massive scalar fields at level four. And they found that the closed string field theory potential 
has a local minimum where both the tachyon and dilaton take positive expectation values. They 

^Several years before the quartic vertex was solved in '7 , Belopolsky [14] managed to calculate the tachyon 
effective potential to order four and he found that it had no minimum. This result, however, didn't take into account 
the dilaton. In fact, it is understood now that the tachyon effective potential doesn't make much sense because one 
cannot integrate out the massless dilaton. 
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also noticed that the depth of this minimum tended to decrease as the level increases, and they 
made the proposition that this vacuum should have a vanishing action density. They supported this 
claim by looking at the low-energy effective action of the tachyon, dilaton and metric. They found 
that if this action has a stable vacuum, its depth must be zero. Although the numerics to level four 
seemed to confirm this claim, a recent computation to level ten [15] shows that, at quartic order, 
the value of the potential at the stable vacuum is actually negative and non-zero. The question 
now is whether higher order terms in the CSFT action can make the shallowness of the potential 
go to zero, or if it stays finite. 

It is clear at this point that we need the quintic term of the CSFT action. In the present 
paper we solve numerically the geometry of the quintic vertex. This is again done with quadratic 
differentials. We present in details the algorithm to solve the Strebel condition and we spend quite 
some time describing the reduced moduli space of spheres with five punctures. It will turn out 
that we can split it into 120 regions, and need to describe only one of them, that we call ^5. We 
undertake the description of the boundary of ^5, and to our pleasant surprise we find that its 
projection on one of its two complex coordinates (corresponding to the two unfixed punctures in 
the uniformizer coordinate) can be described algebraically in terms of an algebraic curve, that will 
be found from the solution of the class of quadratic differentials with two double zeros. The rest of 
the boundary will be solved numerically. After this is done, we can integrate correlators over the 
reduced moduli space. The simplest one is the term with five tachyons. We describe in details how 
we do this integration and how we estimate the uncertainty in the result. 

In order to gain confidence in this result, and a fortiori in the machinery developed and the 
numerical results produced, we must check our algorithm in some way. For this, we compute 
the effective potential of the dilaton to order five. As the dilaton theorem claims, this should be 
identically zero. The term of order five is composed of two terms, namely the contact term that we 
calculate with our algorithm, and the term from Feynman diagrams with vertices of lower order, 
that we calculate with the techniques and results of [11^ [TS] . The cancellation is achieved with an 
accuracy of about 0.1%, falling well within the 0.5% estimated error on the five-dilaton contact 
term. We therefore claim that our algorithm and result for the five-tachyon term, are reliable. The 
computation of terms of higher levels, necessary in order to pursue the study of the nonperturbative 
vacuum of \12\ ITS] , is now possible and will be done in a future publication [16] . 

We think it is of interest to give some orders of magnitude related to the algorithm developed 
here. Its implementation on a computer is done in the C-I-+ language, and the code is more than 
10,000 lines long. The complete and accurate computation of the moduli space and the quadratic 
differentials inside it, takes several days to compute on a desktop computer and generates about 
one GByte of data. Once this is done, the accurate computation of an integral takes several hours. 
Unfortunately, the size of the data makes it impossible for now to express the numerical solution 
in terms of a reasonably short fit, as we did for the quartic vertex [7]. 

The paper is structured as follows: We end this section by a short summary on the CSFT action. 
In Section [2l we construct the quadratic differentials pertinent to the five-point vertex, and we 
study some of their limiting cases in Section [3l We then describe in Section S] how to solve the 
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Strebel condition numerically. The reduced moduli space and the way to compute it are described 
in Section [5j We can then integrate to obtain the first result, namely the five-tachyon contact 
term, this is done in Section O The reliability check with the computation of the dilaton effective 
potentials is done in Section [71 We end with some discussions on the present results and prospects 
in Section [HI 

We now summarize the CSFT action and fix our notation. With a' = 2, the closed string field 
theory action is [U [2] 

S = -^ Q(*|co Qi?|*> + ^ ^> + ^ + ^ + (1-2) 

where Cq = ^(cq — Cc), Qb is the BRST charge and {. . .}, are the multilinear string functions. In 
this paper we won't need to know much about the string field ^, except that when considering 
tachyon condensation we keep only scalars with zero-momentum. Namely 

\^)=t\T)+d\D) + Y,i^i\^i), (1.3) 

i>2 

where ^i, i > 2, are massive scalars, and the first two fields are respectively the tachyon 

|r)=cici|0), (1.4) 

and the dilaton (jl.ip . We will often use the notation V^.^^.^,„^.^ to designate the coefficient of 
il)i^il>i2 ■ ■ ■ tpij^ in the potential (and it is understood that tpi = t and Tp2 = d). 

In the CSFT constructed in [H [11 |3l |3| , the interacting worldsheet id endowed with the following 
geometry: Every external state spans a semi-infinite cylinder of perimeter 27r, and these cylinders 
intersect on a prism. The prisms of a contact term must have the characteristic that all their 
nontrivial closed curves must have length greater than or equal to 2tt, which for the five-point 
vertex is equivalent to all edges having length smaller than or equal to vr. It can be shown that 
all other prisms are obtained from Feynman diagrams with vertices of lower orders. The relevant 
prism for the five-point contact term was first discussed in [3j. It is shown on Figure [H It is made 




Figure 1: The relevant prism of the five-point vertex. The labehng of lengths will be kept as shown. 
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of two opposing triangles, connected with three quadrilaterals. The perimeter conditions on the 
cylinders can be written 



^1 + ^2 + ^3 = 2^ 



4 + ^5 + ^6 = 2^ 



f 7 = ^1 + ^4 - vr 



l9=i3 + h-Tr. 



(1.5) 



Those are five conditions, leaving four independent lengths, matching the real dimensionality of the 
moduli space of spheres with five punctures. All other prisms with five faces would have less than 
four independent lengths and thus correspond to subsets of the moduli space with measure zero. 
Although they don't contribute to the integration over moduli space, some of these prisms will be 
considered in Section [3] because they will be useful to calculate the boundaries of the second kind 
described in Section El 

The usefulness of this particular geometry is that it arises from a quadratic differential ^ that 
has second order poles at the punctures where the external states are inserted, and verifies the 
Strebel condition that its critical graph has measure zero. The ring domains correspond to the 
semi-infinite cylinders and the critical graph corresponds to the prism. 

To be concrete, we will always map the sphere on the complex plane and we will fix the topology 
of the vertex from the beginning, i.e. we will use the labeling of zeros and poles as indicated on the 
right of Figure O Fixing the topology means that we will only consider configurations obtained 
from this one by a continuous transformation and without any two zeros merging. In other words, 
the punctures on quadrilateral faces will always be mapped to z = 0, z = I and z = oo, whereas 
the remaining punctures at and ^2 will always correspond to triangular faces. In these notations, 
the two complex numbers ^1 and ^2 parameterize the moduli space of five-punctured spheres. 

2 Quadratic differentials 

To describe the right geometry and the local coordinates on a punctured sphere, we need a quadratic 
differential transforming like 



under a conformal change of variable. It should be holomorphic everywhere except at the punctures 
zj, where it has poles of order two with "residue" minus one, and has thus the expansion 



If we place a puncture at infinity, as we will, the expansion of (p in the coordinate t = 1/z is 




(2.1) 




-1 



(2.2) 



(2.3) 



It is easily seen that the quadratic differentials obeying these conditions are given by 



0(z) 



z6 + P^{z) + ai z{z -l){z- 6)(^ - 6) + Q2 z\z - l){z - 6)(^ - 6) 

,2 (^_l)2(^_^^)2(^_^2)2 



(2.4) 
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Here P^iz) is a polynomial of order five which is partially determined by the residue conditions 
and will be given by (j2.12p . and oi and 02 do not change the residue conditions at the poles. The 
coefficient of in the numerator must be —1 for the residue at infinity to be —1. The parameters ai 
and a2 will be determined by the Strebel condition, namely the condition that the critical graph of 
(p closes (see [8l[9lll41[7] for more details). It will be enough to know here that the Strebel solution 
gives the quadratic differential needed for the vertex. It can be expressed as the requirement that 
the complex lengths between any two zeros zi and Z2 of the quadratic differential, are real 



Z2 



Im / y/^dz = 0. (2.5) 
J 21 

Solving numerically the Strebel condition will be the object of Section HI For now let us go back 
to P^iz). We will write it in the form 

P,{z) = hz' + biZ^ + P3iz), (2.6) 

where 

Psiz) =b3Z^ + b2Z^ + biz + bo. (2.7) 

The four coefficients of P^iz) can be completely determined, in terms of .^1, ^2, ^5 and 64, by the 
four residue conditions at the finite poles 0, 1, ^1 and ^2- These conditions can be written 

yo ^ ^^3(0) = ^? el 

y, = P3(l) = -l-65-64+(l-6)'(l-6)' 

2/2 = P3{^i) = -Ct-b,fi-b4Ct + ei{^i-l?iCi-^2? 

2/3 = ^'3(6) = -e2-fe5ef-fe4d + d(6-i)'(ei-6)'. (2.8) 

The polynomial P3 that satisfies ()2.8p can be written 

(z-l)(z-6)(^-g2) , z(z-gi)(z-6) ^ 

"^"^'^ = =66 (i-6)(i-6) 

z{z-l)iz-^2) z(z-l)(z-6) .2 

6(6 - i)(6 - 6) 6(6 - i)(6 - 6) ■ ^ ■ ^ 

Now we want to choose 65 and 64 in (j2.6p in such a way that the expression (|2.9p is as simple as 
possible. For this we note, from ()2.8|) . that if 

64 + 656 + 6' = 

and 64 + 656 + ^1 = 0, (2.10) 
then the third and fourth terms in (12. 9p will be simplified. Even better, since the solution of ()2.10p 



IS 



65 = -(6 + 6) 

64 = 66, (2.11) 
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we have that —1 — 65 — 64 = — (1 — ,^i)(l — ,^2)) and thus we see, from (j2.8p . that the second term 
wiU be simphfied too. And obviously, since yo = Ci'^ii ^^^o have that the first term is simplified. 
Thus, making the choice (j2.11|) for 65 and 64, we have that 



P^{z) = -sz^ + tz^ + (-S + v{s - 1)) + {t + s^ + v{l-s- t)) z^ + t{v -2s) z + t^ , (2.12) 



where 

s = 6 + 6 
t = 6 6 

V = (6-6)' = s'-4t (2.13) 
are symmetric expressions in 6 and 6- 

We note, for future reference, that the derivatives of (/)(z) with respect to oi and 02 are 

d(l){z) _ -1 

dai 2(2 - 1) (z - 6) (^ - 6) ' 

d(j){z) _ -1 

~ (z-i)(z-6)(^-6) ■ 

The most uniform prism 



(2.14) 



We cah the most uniform prism, the prism (Figure [T]) which has two equilateral triangles. The 
lengths of the edges of the triangles must therefore be and the three edges connecting the two 
triangles must have lengths ^. We want to describe the quadratic differential of this configuration 
in the z-plane, where we will put three punctures at respectively zero, one and infinity. As already 
mentioned, we have to start by deciding which faces will correspond to these three punctures. Since 
we are later going to consider the subgroup of PSL(2, C) conformal transformations that permute 
the aforementioned punctures, it will be simplest to map the three quadrilateral faces on z = 0, 
z = 1 and z = 00. We must then stick to this choice because it would be very difficult for the 
numerical algorithm to switch between configurations of different topologies; this will be discussed 
more extensively in Section [5l For the purpose of finding the quadratic differential of the most 
uniform prism, it is easier to work in another coordinate w, where the symmetry of the prism is 
more visible (see Figure[2]). In the w coordinate, (jy^w) has three poles at the vertices of an equilateral 
triangle, respectively w = —1, w = e's" and w = e~~^ , which correspond to the quadrilateral faces. 
One pole is at its center w = 0, and the last pole is at infinity. By contemplating the left half of 
Figure [21 we can immediately write the ansatz for the zeros Wi, i = 1, . . . , 6, of (j){w) 

wi = P , W2 = e 3 (j , W3 = e 3 p 

1^4 = 7, W5 = e 3 J , wq = e 3 •y , (2-15) 
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Figure 2: The critical graph of the quadratic differential corresponding to the most uniform prism, in the to- plane 
where the symmetry is obvious, and in the z-plane in which the three quadrilateral punctures are at the standard 
points 0, 1, oo. 



where /3 and 7 are positive real numbers and 7 > /3. The quadratic differential for this configuration 
is thus 

, , _ —{w — f3){w — j){w — e^f3){w — e~^j){w — e~^P){w — e~^j) 
cpyW) - — 



{w + ly {w — (^ui — e 

{w^ - P^) {w^ - 7^) 



^2 (^^3 _j_ 

The residue condition at w = —1 and w = are respectively 



(2.16) 



(/?3 + l)(7»ti)=l , p^y^^i. (2.17) 
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The solution to the system (I2.17P with the constraint 7 > /? > is 



1 , _ 1 

3 



,^{i^y. .-(i±^) . (2.1S) 

Now we express the quadratic differential in the z coordinate. We want to map the three 
quadrilateral punctures to z = 0, z = 1 and z = 00. Namely z = h{w), where 

h{-l)=0 , h{ef) = l , /i(e-f) = oo, (2.19) 

whence 

z = h{w) = e- (2.20) 

w — e 3 
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which has the inverse 



Remembering that 



we find 



w = h-\z) = '~''\^' . (2.21) 
z — e 3 

(dw \ ^ 
^ ) ' ^^-^^^ 

, . _ - + - + - 3z + 1 



where the poles and .^2 are 

^1 = h{oo) =63 

6 = /i(0) = e-"^ . (2.24) 

Now we can easily determine the parameters ai and 02 by comparing ()2.4p and (j2.12p with (j2.23p . 
We find 

ai = a2 = -2. (2.25) 

At last, the zeros of (piz) are given by 

Zi = h{wi) , i = l,...,6, (2.26) 

and their positions are shown on Figure [21 It is important to know an exact quadratic differential 
for the numerical algorithm to start with. It will then cover all the reduced moduli space by 
successive deformations of this solution, each deformation being relatively small for the Newton 
method to converge (more on this in Section . 

3 Limits of quadratic differentials 

One of the most complicated problems in the computation undertaken in this paper, is to describe 
the boundary of the reduced moduli space. As was already noted in [3j, there are two kinds of 
boundaries. When the length of one of the triangle edges is vr, we are on a boundary of the first 
kind, corresponding to the situation in which the vertex can be built as a Feynman diagram with 
a propagator of zero length. There are also boundaries of the second kind, when one of the lengths 
£7, is or £9 is zero. These are internal boundaries, they connect prisms with different assignations 
of punctures to faces. They are not boundaries of the whole reduced moduli space, but they are 
boundaries of the region A5 that we will consider by keeping a fixed assignment of punctures to 
faces (see Section [5]). We will also need to consider intersections of boundaries, in particular when 
two edges have vanishing lengths. The description of the moduli space will be done in Section \5\ 
The problem for now is that when a length vanishes, the quadratic differential has a double zero 
and it becomes numerically ill-conditioned. The goal of this section is precisely to deal with these 
limit cases. We can get rid of the numerical difficulties by some analytical work. Actually in the 
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case of two double zeros, the Strebel differentials can be explicitly described in terms of an algebraic 
curve, reducing the numerical work to finding the roots of a polynomial of order six§| 



3.1 Quadratic differentials with one double zero 
The general case 

When we compute the boundary of the reduced moduli space, we are led to consider configurations 
where the quadratic differential has a double zero, or in other words when an edge collapses to zero 
length. We have to distinguish whether the collapsed edge is an edge of a triangle or one of the 
three edges connecting the two triangles. In the first case we end up with a face with only two 
edges, and for this configuration to be in the reduced moduli space, they need to have length vr. 
The remaining six lengths are constrained by four residue conditions, and we are therefore left with 
two real degrees of freedom. In the second case, we have eight lengths constrained by five residue 
conditions, which leave us three degrees of freedom. We are thus going to consider this case only as 
this will describe some boundaries of the reduced moduli space. Once we have fixed the topology 
and the labeling of zeros and edges (Figures [U and [2] and Equ. lll.Sp ) in the z coordinate, we must 
describe separately the cases where £7, or ^9 respectively, are zero. Since we have now only three 
real degrees of freedom we can fix ,^1 and only one real component of .^2) for example 



where C is a given complex number of unit norm. Once these quantities are fixed, the quadratic dif- 
ferentials with one double zero, satisfying the residue conditions have three real degrees of freedom, 
which are for example, Ren, Imn and (iCI?2)) where u is the position of the double zero. 

We will thus write 




(3.1) 




(3.2) 



where the cubic polynomial 




(3.3) 



will be completely determined by the following residue conditions at the finite poles 




(3.4) 



^ Other analytically solvable limits of quadratic diflterentials with four and five poles have been studied in the 
context of open-closed duality in [17] . 
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The polynomial P3{z) that satisfies p.4p . can be written as in ()2.9p . Now if we define 

2/0 



^^0 

Vi 

V2 



-66 

yi 



(i-ei)(i-6) 

y2 

6(6-i)(ei-6) 



= 6(6-1X6-6)' ^'-'^ 



we have that 



^3(2) = {VQ + VI + V2 + V3) - {{1 + S)V0 + SVl + {I + C2)V2 + {1 + CM z"^ + 

+ {{s+t)vo + tVi+C2V2+ClV3) Z+ (-] , (3.6) 



where again, s = 6 + 6 t = 66- And the quadratic differential is completely determined 
once we give ^1, 6 and u. 

The regular pyramid 

As we did for the regular configurations, we want to calculate explicitly the Strebel quadratic 
differentials with one double zero in the most symmetric case, a pyramid with a square base (with 
edges of length |) and four triangles with edges of lengths ^, ^ and Again, it is easier to 
solve the quadratic differential in another coordinate w, where the symmetry is obvious, and then 
map it to the z-plane. We will have three different mappings hi, /i2 and /13, whether the vanishing 
length is ^gj ^9 oi" ^7 respectively. In the w coordinate we set the double zero u = and four poles 
around it, namely at 1, i, —1 and —i. The last pole, corresponding to the base of the pyramid, is 
at u) = 00. By symmetry (see top of Figure ED, we can then immediately read off the ansatz for 
the zeros Wj, i = 1 . . . 4 of 4'{w), namely 

wi = ae i , W2 = ote 4 ^ = ae 4, , = ae 4, , (3.7) 

where a is a real positive number. We have thus 

{w^ - 1) 



^iw)=- (3.8) 



The residue condition at the pole w = 1 gives us 

1 + /si /X 

__ = 1 ^ a = (15)4. (3.9) 

We now map it in the z-plane to the configuration (lower left of Figure [3]) which has = 0. The 
map is 

z = hi{w) = '^^. (3.10) 
w + 1 

This maps the double zero to 

mW = /ii(0) = -1, (3.11) 
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Figure 3: The regular pyramid in the lu-plane and its three mappings to the z-plane with respectively = 0, = 
and ^^ = 0. The double zero is marked with a small circle. 



and the poles are 

cf^ = /ii(i) = i , e5^) = /ii(-i) = -i. (3.12) 

And the quadratic differential is 

For the configuration with ^9 = (lower middle of Figure [3|) , we have 

^2{w) = \^{W + 1) 
(2) 1 

2 

cP = ^(i + o , = ^(1-) 

</.(2) z = ^ 1 ? ^ . 3.14 
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And for the configuration with ^7 = (lower right of Figure [3]) , we have 



2 



u; + 1 



u 



(3) 



= 2 



(3) 



1 + i 



= l-z 



1 




</)(3)(z) 



z^z - 1)"^ {z - {1 + i)f (z - {1 - i)) 



(3.15) 



It is very useful to have these three exact quadratic differentials to start the Newton method 



when plotting the boundaries of the second kind. 

3.2 Quadratic differentials with two double zeros 

As we will see in Section [5l it is important to be able to construct quadratic differentials with two 
double zeros. In the present section we will see that we can actually solve the Strebel condition 
in terms of an algebraic curve, i.e. everything can be done algebraically except for a root of a 
polynomial of order six, which must be found numerically. 

From Figure [1] and the lengths conditions (jl.Sp . we see that when the quadratic differential 
has two double zero, one face must have only two sides (which must therefore have length it). It 
can then be further deduced that only one length is free, with value i. The other lengths are 
respectively £, n — i (two edges) and ir. We now look at the quadratic differential in the w-plane 
where the puncture at infinity is attached to the face with two edges, and the punctures of the 
other two faces which have these edges as a side are mapped to w = —1 and w = 1. It is readily 
seen that in this coordinate, the critical graph must have the symmetry w —w. This is shown 
on Figure HI The two remaining punctures will then be at, say, w = a and w = —a. The two 



Figure 4: The configuration with two double zeros in the «;-plane with the particular value £ = 0.8. The length of 
each edge is shown. 

double zeros are at w = u and w = —u, and the two simple zeros are at w = b and w = —b. Given 
the length £, we want to determine the three complex numbers a, u and b. 




I w 



TT 



TT 
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The quadratic differential Lp = 4>{w){dw)'^ for this configuration is 



ct){w) 



{w — u)"^ {w + u)^ {w — b) {w + b) 



(3.16) 



(w — 1)2 (w + 1)2 (tu — a)2 (w + ay {w'^ — 1)^ (t/;2 — a'^y 

The residue condition at infinity is automatically satisfied by the expression (j3.16p . and by sym- 
metry the residue condition at w = — 1 is the same as for tt; = 1, and the condition at w = —a is 
the same as the one at w = a. We thus have only two independent residue conditions, at w = 1 
and w = a respectively, which read 

(l-n2)2(l-62) _ 



4(l-a2)2 



[a 



1 



(3.17) 
(3.18) 



4a2(l -q2)2 

We are now going to solve the Strebel condition. For this, we need to write that the complex 
lengths between any two zeros are real. By symmetry, this is automatically satisfied if one length, 
say the length between b and n, is real. We therefore impose 



where 



£{b,u) 



" (w;2 



tx2) VP 



(1-^1-2 



dw . 



(3.19) 



(3.20) 



And in total we have two complex equations and one real equation to determine six real parameters, 
we thus have one free real parameter for the configurations with two double zeros. We take this 
parameter to be £ G [0, 7r/2]. The other values £ G [7r/2,7r] are trivially related to the first case by 
complex conjugation. 

The integral in (|3.2U|) can be calculated by doing the substitution y = \ 1 



fe2 



(y2 - 1) (y2 _ (1 _ 52)) (^y2 _ ( 1 _ ^ 

dy 




whereby 
dy 



(3.21) 



To go from the first to the second line we have made use of the residue conditions (j3.17p and (j3.18p . 
So the length between b and u is 



£ = - + - In 
2 2 



1 



62 

1/2 



1 



62 



62 
1^2 



4 + 1 



62 

7/2 



fe2 
1^2 
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It is natural to make the definitions 



s^\ll-^ , t^V^ , v^^Jl-^. (3.22) 



We thus have 



{s + l){s-t)^s + vy ■ ^ ' 
We note that the equations (j3.17p and (jS.lSp give us a simple relation between t and v 

t-v = -hv. (3.24) 

It will be convenient to define 

T = tv. (3.25) 
The following identities follow directly from (|3.24p 



We also have from (|3T7D . (|3:22D and (15:26]) 

2 _ 3r(r - 4) 

* ~ ~ 4(1 + T) ' 

and 



(3.27) 



ib ih 

u = ^== , « = ^^- (3-28) 
vs — 1 vf^ — 1 

Let us now rewrite (j3.23p as a polynomial equation in s, whose coefficients are expressed in terms 

of T and i with the help of (fOHIl 

- i cot(^) - (1 + T)/ - i cot(£) - + (^T^ + 2r^ - r2 = . (3.29) 

Using (j3.27p this becomes 

,, = _lu„(f)l!E±il(£zl) (3.30) 

3 ^ ^ (T - 4)2 ^ ' 

Squaring and using again (I3.27p . we finally get a polynomial equation for T. 

Pt{T) = -coi^{l)T{T-Af + A^{T + lf(T-^^ =0. (3.31) 



9 , 

We immediately see the solutions of this equation for the particular values £ = 7r/2 and I ^ Q. 
Namely T = 16/9 and T = 4 respectively. For generic this equation must be solved numerically, 
following the branch T{1 = 7r/2) = 16/9. 
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We would like to characterize precisely the branch of the solution of (j3.3ip that we must choose. 
We start by showing that the branches cross only at the points T = ^ and T = 4. Branches cross 
at multiple zeros, when Pe{T) and P^iT) vanish simultaneously. We have 



p;(r) = (3r - 2) ( -2 cot^{£){T - 4)^ + - 48(r + 1)^ 
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We define 



A = - cot\i){T -4f 
so that we have the system 



B = 48(r + 1)^ ( T - y 
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= Pg(T)=T{T-4:)A + {T+l)\^T- — ] B 
= P!,{T) = {3T -2) (^2A+^B^ . 



Let us look at the second equation. If T = g, (|3.3ip tells us that cot (^) 
have B = —^A, and the first equation tells us that either T 



forT 



(3.32) 
(3.33) 

(3.34) 
(3.35) 



-1. If T 7^ |, we must 
-16. In both of these 



cases we see from (I3.3ip . that cot'^(^) = — 1. We have thus shown that the branches do not cross 
when i € (0, ^). From this fact we can completely derive the topology of the branch diagram (by 
which we mean the locus in the T-plane formed by all roots of Pi{T) for all £ S [0, §]). Since the 
polynomial Pe(T) is real, the branch diagram must be symmetric under complex conjugation. Five 
branches must start (at £ = 0) from T = 4 and one must start from T = 0, whereas three branches 
end (when £ = ^) at T = —1 and two must end at T = Also when £ ^ ^, cot'^{£) tends to zero 
and we have a solution 

48 

T ^ -7TT^ +00 , 



cot'{£) 

and therefore one branch must end at +oo. It is then clear that the only topology compatible with 
all these observations is the one shown on Figure O The branch corresponding to the conventions 
of Figure H] is the one drawn with a bold line. It can be characterized by the fact that it belongs 
to the rectangle ^ < ReT < 4 and —0.6 < ImT < 0. We will thus define the function T{£) by 



T{£) is the unique solution of < 



cot2(^) T{T - 4)5 + 48(T + 1)^ 
16 



9 



9 



< ReT < 4 



-0.6 < ImT < 



(3.36) 



From (j3.22p . (j3.26p and ()3.28p . we can write the expressions of b(£), u{£) and a{£) in terms of r(^). 
The only one we will need is a{£) 



a{£) = i 



2T{£) - 1 - |r(^)2 



(3.37) 
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Figure 5: The branch structure of Eau. (|3.3ip . The branch chosen is the bold one. 



This expression will be useful later do describe the boundary of the projection of the reduced moduli 
space on the ^i-plane. We already note the special values 

27 

«(0)=2-V5 . „(V2) = ^^^^_^i. (3.38) 
4 Solving a quadratic differential numerically 

We now want to describe the numerical algorithm to solve a quadratic differential. We will focus on 
the regular quadratic differentials, but the techniques can be applied, with only trivial modifications, 
to quadratic differentials with one double zero. We will then recall how to compute the mapping 
radii, and we'll shortly discuss how we compute derivatives with respect to and ^j. 

The central problem is to find, for given poles and the quadratic differential (i.e. its 
parameters ai and 02) that satisfies the Strebel condition. In other words its critical graph must 
close, with its zeros being linked with horizontal trajectories. We remind that horizontal trajectories 
are defined by the condition (j){z){dz)'^ > 0, so the Strebel condition is equivalent to the condition 
that all the complex lengths between zeros, Zi and Zf, are real. 



lme{zi,Zf)=lm y^dz = , i,/ = l,...,6. (4.1) 

J Zi 

We have nine lengths (see Figure [2|) constrained by the five residue conditions; this leaves us four 
independent lengths, which we can choose as iQ, £4, £g and £2 (as labeled on Figure [2]). We thus 
have to solve a system of four real equations of four real unknowns 

f(x) = 0, (4.2) 
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with f = {lm£Q,lm.£4,lm.£g,lmi2)^ and x = (Re ai, Im ai, Re 02, Im 02)"^. We use Newton's 
method 

^^+1 = - h:; • f (^^) ' (4-3) 



which converges when the initial guess xq is not too far from the solution. The Jacobian 5f /5x 
can be expressed in terms of derivatives with respect to Oj 



5f 

dx. 



( Im5ai4 Re5ai4 Im^aa^ Re^aa^ \ 

\m.da^l/4 Re5a^£4 Im5a2^4 Re 502^4 

Im^ai^g Re^ai^g ImSaj^g Re Sag 4 

V Im5ai^2 Re5ai^2 \m.da^i2 Re(9a2^2 / 



(4.4) 



We continue this section by explaining how to accurately compute numerically the complex 
lengths and their derivatives with respect to ai. And we will continue by describing the computation 
of the mapping radii, and of the derivatives of with respect to and ^j. 

4.1 The complex lengths 

The quintic numerical computation turns out to be very much harder than the quartic. In the 
quartic computation [7], there were two weaknesses which had negligible consequences. First it 
was hard to tell if the path of integration (which was chosen to be a straight line) was going on 
the right side of the poles, hence a possible ambiguity of 27rn in the length. Combined with the 
sign ambiguity of the square root, there was a potential problem when the length was nearly vr. 
Although we could go around this ambiguity in the quartic case, for the quintic calculation it would 
be catastrophic. Second, we used the coded "continuous square root", which remembers its last 
evaluation and tries to detect if the branch cut has been crossed between the last two evaluations. 
Again this is not good enough in the quintic case; we would inevitably encounter situations when 
the continuous square root fails to detect a branch cut, again with catastrophic consequences. 

To solve the first problem, we make a conjecture based on observation that the integration path 
between two zeros Zi and z^ along the critical trajectory C that connects them, can be continuously 
deformed to the path shown on Figure [U to "Pj to without crossing any pole (other than pj). 
Here is a pole such that the zeros Zi and Zj are on the boundary of its ring domain. 



l{zi,z^)= \ ^/^dz= / 70(1) dz+ / ^/^dz+ / ^/^dz. 

Jc Jzi Js Jz2 

For numerical reasons that will become clear in a while, we split the integrals along the straight 
lines in two, by setting mi = ^{zi + pj) and m2 = ^{zf +Pj), and writing 

mi t-zi /• /•m2 rzf 



£{zi,Zf)= / ^/^dz+ / ^/^dz+ / ^/^dz+ / ^/^dz+ / ^/^dz. (4.5) 

J Zi J mi Js J Z2 J m2 

Now we subtract the poles at z = pj in the integrals from to z^, and put them back together 
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Figure 6: The critical trajectory C can always be deformed continuously to the contour integration shown with 
solid lines. 



with the integral over S", and we take the hmit zi,Z2 — Pj ■ 



mi /•m2 



izi,Zf) = / ^y (p{z) dz - / ^J^{z)dz 



+ r fv^- dz - r fy^R - ) dz 

+ lim ( / \/ 4>{z) dz + r [ — ^ — dz — r [ — ^ — dz] 



(4.6) 



ki-Pjl=|2^2-Pjl 



where r = ±i is the residue of \/ ^{z) sX z = pj. It is now time to define y^0(i). We want it to be 
continuous on the integration path, in other words, we don't want to cross any branch cut. The 
idea is to calculate the intersection of the integration path with the branch cut of the conventional 
square root, and then change the sign according to the number n{z) of branch points that were 
crossed. Since ^(z) is given by an expression of the form 

= (4.T) 

we write 

^/^-^^^3(_l)n(.), (4.8) 

where the square root on the right hand side is the conventional square root with branch cut the 
negative real axis. We ignore r] for now and focus on n{z). When zq is on the integration path, 
say [zi,pj], we define n{zo) to be the number of times that the conventional square root ^/—ri'^P{y) 
crosses its branch cut when y G [zo,^^], In particular n{pj) = 0. We start by parameterizing the 
path with t e [-1,1] 

z{t) = ^^t + ?i±^^az + b, (4.9) 
n{zo) is then given by the number of solutions U of the equation 

lm{-r,^P{z{U)))=0, (4.10) 
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satisfying the conditions 

Re (-??2p(z(ti))) < and to < U < 1 , (4.11) 

where z^Iq) = zq. It is then clear that (j4.8p will be continuous. Let us now look at the complex 
number rj in (j4.8|) . We want to choose it in such a way that the residue of y(j){z) at z = pj , which 
is ibz, takes the negative sign. Remembering that the "residues" of 0(z) at the poles are —1, we 
see that P{pj) = Yli^^jiPj — PiY — {^'{Pj)Y- We thus require that -q satisfies 



(4.12) 



riD'{pj) 
which is the case if 

lra{rjD'{pj)) >^ or {lxn{r]D' [pj)) = {} and Re (7?L>'(pj)) < O) . (4.13) 
We will therefore take (demanding also that rj has unit norm) 

' D'{pj) ' ^^-^^^ 

where Q € (— ^, ^) can be chosen arbitrarily, as long as as it is not zero. Indeed if = 0, we would 
have —T]^P{j>j) > 0, and (j4.10p would have the solution t = 1. This must be avoided because, 
due to numerical uncertainties, we would find a solution t = 1 it e and it would be hard to decide 
whether to count it in n{z) or not. li 6 ^ 0, we avoid this problem. 

We can now continue the computation (j4.6p by setting r = —i. In the limit considered, the 
integral over the arc of circle S is just a bit of the residue, namely 



\/ 4>{z)dz = ai{—i) = a , (4-15) 

where a is the angle at the tip pj of the triangle ZiPjZf. We thus have 

^/^dz = arg* ~ , (4.16) 
s Zi- Pj 

where arg* is the argument function with range [0, 27r), in other words, its branch cut is on the 
positive real axis. This is necessary because a might be greater than tt, when pj is on the right of 
the straight line from Zj to z/ but still on the left of the path C. The last two integrals in ()4.6p give 

lim ( —i I dz + i I dz 



|zi-pl = |z2-Pjl 



lim log I — — — ^ ] = — ilo, 



Pj 



■Pj 



(4.17) 



where in the last equality, we used the fact that zi is on the segment from Zi to pj and Z2 is on the 
segment from Zf to pj, and the argument of the log is therefore a positive real number. Combining 
([H^ and (fiTT]) . we see that the last line of (fOj) is 

-^log*^^, 

Zi - Pj 
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where log* z = log 1^1 + i arg* z. In total we have therefore 

i{zi,Zf) = / yQ^dz+ — \ dz 

J Zi J mi \ ^ Pj / 



'■m2 fPj 



Zf Jni2 



y^dz+ r dz] -Hog* (4.18) 



z 



Pj J Zi- Pj 



We shall now explain why we have split in two the integrals along the straight lines. For this we 
need to make a parenthesis into the Gaussian quadrature formulas (see for example |18j ) which we 
use to compute numerically the integrals in (j4.18p . Those are very useful when integrating, over 
the finite interval [—1, 1], a function g{t) that behaves like (1 — t)" near t = 1, and like (1 + t)^ 
near t = —1, where a,f3 > —1. If a or /3 are not integer the function g{t) will have singular 
derivatives, and a quadrature formula like the trapezoidal method or the Gauss method with no 
weight (which is adapted when a = (3 = 0) will give inaccurate results. In that case, we need to 
use the Gauss-Jacobi quadrature formula 

1 N 

/ {l-tni + tff{t)dt = y^H,fia,)+E, (4.19) 

where 

g{t) = {l-tni + t)f'f{t), (4.20) 

and f{t) is infinitely differentiable with all its derivatives bounded on [—1, 1]. The abscissas aj are 
the roots of the Jacobi polynomial JAr(t;a, /?), the weights Hj are given by 

^ 2A^ + a + /3 + 2 r(iV + a + l)r(iV + /? + l) 2^ 

and the error E is 

^ ^ TiN + a + l)T{N + (3 + l)r{N + a + P + l) iV!2^^+"+/^+^ ^ 

(2iV + a + /? + 1) (r(2iV + a + /3 + 1))2 (2iV)! > I > J- 

(4.22) 

We take = 15 for all our integrations, which is enough to obtain a typical accuracy of fourteen 
significant digits. 

Now suppose that we didn't split the integrals in (14. 5p and tried to integrate 



At t = —1, the first term behave like {t + 1)^^'^, but because of the second term, we cannot write 
the integrand in the form ()4.20p . This is why we had to split the integrals in (j4.5p . In this way, the 
first and third integrals in (|4.18|) are of the form (|4.20|) with a = and /3 = i, while the second 
and fourth integrals have a = (3 = 0. 

Note that there is an alternative to the splitting, which was used in [7J. It is to subtract to 
(/>(z) the expression 



-i y/TTT 



V2 {z{t)-pj) 



22 



instead of —i/{z{t) —pj)- In that way we can factor out the square root along the whole path. But 
we find the splitting method easier, because the expression 



Z-pj 

is simpler when one simplifies it by explicitly canceling the poles. What we mean is that this 
expression is prone to numerical errors, because when z is close to pj we are subtracting two large 
and almost equal numbers, a cancellation error. So we need to manipulate this expression. Namely 

I i ^ V-r]^P{z) ^ i _ ^J-r]^P{z)+ir]Dj{z) 



Pj rjD{z) z ~ Pj rjD{z) 

^D^{zl^^ (^.■(^) {V^l^)-^nD,{z))y' , (4.23) 



where we have defined 



D{z 
z-pj 

We can cancel the apparent pole in the factor 



Q(.) . 2MLzIMl ,4,25) 

z-pj 

because the residue condition of the quadratic differential at its poles implies that P{pj) = Dj{pj)^, 
so the numerator of (j4.25p is a polynomial with a factor z — Pj- Thus Q{z) is a polynomial, and 
the right hand side of ()4.23p is well defined, analytically as well as numerically. 

In order to use the Newton method, we still need to calculate the derivatives daji{zi, Zf). Since 
(j){z) vanishes at the endpoints of the integration path, we can simply differentiate inside the integral 

I-.,., . (4.20) 

With the expressions (|2.14p for the derivatives and the definition (|4.8p of the square root, and 
deforming the integration contour as before we find 

^£{z„Zf) = --{ ^ dz- ' dz] 

-^l{z,,zf) = / ; ' dz- ^ ) dz) . (4.27) 



These integrals can be evaluated accurately with a Gauss- Jacobi formula with a = and /5 = — i 



2' 
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4.2 The mapping radii 

We recall that the mapping radii pj appear as a multiplicative factor in front of the local coordinate 
Wj in the power expansion of the maps hj from the local coordinates to the uniformizer z (or 
t = 1/z for the puncture at infinity) (see Equ. (j7.26p l. While all other terms in the expansion can 
be determined by expanding the quadratic differential, the mapping radii can't. Instead, they must 
be computed by integrating sj (t>{z) from a point on the boundary of the ring domain of the pole, 
to the pole pj , and subtracting the singular logarithm (see [21 E] ) 



log Pj = hrn ( Im / 4>{z)dz + log 



a 



Zi - Pj 



(4.28) 



\Zi -Pj\ 

and we take Zi to be one of the zeros on the boundary of the ring domain. The square root is again 



defined as in (j4.8|) with r] as in (j4.14|) . so the residue of (p{z) at z = pj is —i. This integral is very 
similar to that for the complex lengths, so we know how to evaluate it numerically. First we break 
the integration in two and regularize the integrand near the pole 



log Pj = lim I Im / (f){z)dz + Im 



Pj+ea 



zl 



dz + log e — Re / 



1 



Z-Pj 



-dz 



where for definiteness we take zi = \{zi +Pj)- Now we can easily get rid of the limit 



log Pj 



lim I Im 

Zl 



Zl 



\l 4>{z)dz + Im 



Pr 



VHz) + 



Im / (p(z)dz + Im 



Pi 



Zl 



VHz) + 



Pj 



Z-pj 

dz + log 



dz — log 
Zj - Pj 



ea 



Zl - Pj 



+ loge 

(4.29) 



The integrand of the second term is manipulated as in ()4.23p - (I4.25p . and both integrals can be 
evaluated accurately with Gauss-Jacobi quadrature formulas. 



The pole at infinity must be treated separately; we must work in the coordinate t = 1/z. The 
radius is then 



log = lim Im / \/ (j){t)dt + log e 



l/z^ 



(4.30) 



where Zi is a zero on the boundary of the ring domain at infinity. One could of course rewrite the 
integral in the z coordinate 



\/ (p{z)dz , 



but we prefer to integrate numerically over a finite interval, so we will stay in the t coordinate. We 
have thus 



log P5 



lini ( Im / ' \/(t){t)dt + Im / f vw) + -j dt + loge - Re I jdt 



1 



2zi 



Im / " ^/W)dt + Im / ( yW) + -]dt + log 



2z, 



(4.31) 
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The integrand of the second term can be manipulated again, as in (j4.23p - (j4.25p . in order to get 
rid of the apparent singularity at t = 0. It is understood that before doing so one has to express 
(j){t) in the coordinate t. 

cP{t) = cPiz = t-')^. (4.32) 

4.3 The derivatives 

In order to compute amplitudes involving not only tachyons, we need the derivatives and 
We don't have expressions for these quantities that we could directly evaluate, so we need to 
estimate these derivatives. We are using a Richardson's formula of order four, which for a function 
/ differentiable at least five times, reads 

fix) = ^ (fix -h)- 8/(x - h/2) + 8/(x + h/2) - fix + h)) + . (4.33) 

As we are working with numbers in double precision (which is usually around 16 significant digits) 
we see from this formula that the most efficient h is approximately h = 10~^, which should give 
a result accurate to about twelve significant digits; taking a smaller h would reduce the precision 
because of cancellation errors. If we write = Xj + iyj, we have d^^ = ^idx^ — idy^) and 5|. = 
^idxi +idy.), for each Oj we thus need to evaluate the four derivatives and j = 1,2. Each 
derivative calculated with (14.330 requires four evaluation of the function, so in total we need to 
solve the Strebel differential at sixteen points in order to evaluate accurately its derivatives at a 
given point. Note that for the quartic vertex [7], the derivatives are calculated ([HI [I2l [15]) by 
differentiating the fit a^^iCC)^ but for the quintic we haven't been able to write a reasonably short 
fit. 

5 The reduced moduli space 

In this section we describe the boundary of the reduced moduli space Vo,5 (which we sometimes 
call simply moduli space, without ambiguity since we always consider the reduced moduli space), 
and the methods to evaluate it numerically. To represent this four-dimensional moduli space, we 
first study its projection on the ^i-plane, this can be done algebraically. Then at every point of the 
projection, we have to describe the two-dimensional section in the ^2-pla'ne; we do this numerically. 

We naturally want to reduce the numerical part of the problem to the minimum. For this, we 
need to figure out what is the smallest part of Vo,5 that we need to describe numerically and still be 
able to integrate over the whole Vo,5. We can do two successive partitioning of the moduli space. 
Firstly, let five string states \^i), i = 1, ... ,5, scattering on a five-punctured sphere, given by a 
particular point m of Vo,5. We can draw, on the sphere, the critical graph of the Jenkins-Strebel 
quadratic differential uniquely determined by m. For almost all m (i.e. all m E Vo,5 except for 
the ones belonging to a subset of measure zero), the graph will delimit a prism consisting of two 
opposing triangles and three quadrilaterals (see Figure [I]) ; we can thus partition the moduli space 
into ten regions, according to which three of the five states belong to the quadrilateral faces of 
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the prism. We conformally map the sphere onto a sphere, requiring that the three quadrilateral 
vertices are mapped to the standard points z = 0, z = 1 and z = oo (there are six such maps, 
but which one we choose is irrelevant as will become clear below). Secondly, as in the case of the 
quartic vertex OE], we can use the six PSL(2,C) transformations that permute the three points 
{0, l,oo}, and complex conjugation, to partition further the ten regions of Vo,5 into twelve parts. 
In total, the moduli space is thus composed of 120 parts, and we need to describe only one of them, 
which we denote ^5. The integration over the whole reduced moduli space can then be written 

dAi A . . . A dXi{^\b{vx,) . . . 6K)l^l)l*2)|*3)|^'4)|^5) = (5.1) 

Vo,5 

F(^'i,^'2,^'3|^'4,^'5) + permutations^ + c.c. , 

where Aj are real coordinates of the reduced moduli space Vo,5, and the surface state (S| corresponds 
to the 5-punctured sphere given by the parameters Aj. The antighost insertions b{vx-) are not needed 
in this section, they will be given in Section 17. 2i We denote 

F(^i, ^2, ^'3|^'4, ^-5) = dAi A . . . A dX^mbivx,) . . . 6(t;A4)l^i)l^2)|^4)|^5)|^3) , (5.2) 

and it is understood that the first three states are inserted on the quadrilateral faces (i.e. the 
punctures zi = 0, ^2 = 1 and = 00). The integrand on the second line of (j5.1|) is the sum of 
the ten permutations of F{^i, ^'2, ^'3|^'4, ^'5), each permutation being determined by which states 
are assigned to the first three arguments (i.e. the quadrilaterals faces). We are then left with six 
integrals that can all be written as integrals over ^5 after pulling back their integrands through 
the aforementioned PSL(2, C) maps, which are explicitly 

1 1 1 z 

z ^ z , z ^ - , z ^ 1 — z , z ^ , z — > 1 , z ^ . (5.3) 

z 1 — z z z — 1 

The six other integrals, over the complex conjugates of the six domains, can be trivially pulled 
back to integrals over the six domains because the transformation is simply complex conjugation. 
This contribution is therefore the complex conjugate in (jS.ip . 



Now we need to describe ^5. Let us define Vq 5'^'°°^ to be the part of Vo,5 whose three quadrilat- 
erals correspond to the punctures z = 0, z = 1 and z = 00. Its projection on the ^i-plane is shown 
on Figure [71 The six mappings ()5.3p . and complex conjugation, transform simultaneously ^1 and 
^2, we could thus choose to fix either one in a given region. We choose to keep ^1 in the region 
A^^\ characterized by being above the real axis, on the left of the line Re^i = | and outside the 
unit circle centered on zero (see Figure [7]); and ^2 is unconstrained by these mappings. So is 
the projection of A5 on the ^i-plane. 

We will now describe the remaining boundaries of A^\ namely the two curves Bi and 82- In 
order to do this we must anticipate a little bit on the computation of the two-dimensional sections 
of ^5 (the regions in the .^2-plane with fixed ^1). To make things clear, if we parameterize the 
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Figure 7: The projection of Vq 5'^'°°'^ on the ^i-plane. The darker region is the projection of As- 

four-dimensional moduli space by the four real coordinates (Re ^i,Im^i, Re ^2,1111^2)) the section 
is the subset of the ^2-plane for which (Re^i,Im^i,Re^2,Ini^2) is inside the moduli space. We 
present now on Figure O a picture of ^5, i.e. at various points ^1 of ^5^^ we show the section 5^^. 
The boundaries of these sections are of two kinds, the first kind is when an edge of the prism has 
length vr, while the boundaries of the second kind are found when an edge has vanishing length 
(these are drawn with a thick line on Figure [8]) . It is clear that when ^1 approaches the boundaries 
Bi or B2, the section S^-^ must shrink to a point. As we will see below, the sections near Bi have 
two boundaries of the second kind, so when this is reduced to a point, we must have simultaneously 
two edges with length zero. Near B2, the sections have only one boundary of the second kind, but 
three boundaries of the first kind. So on B2 we must have simultaneously one length zero and three 
lengths vr. As can be checked on Figure [1] and Equ. (|1.5p by exhaustion of all possibilities, this 
implies that two edges have vanishing lengths. In other words, a sphere on the boundary Bi or B2 
is described by a quadratic differential with two double zeros. We now understand the importance 
of the quadratic differentials discussed in Section 13. 2^ and we are glad that they could be solved in 
terms of the algebraic curve T{ll.) Equ. (|3.36p . 

What we must do now is to list all the conformal and anti-conformal maps z = h{w) such that 
the poles {—1, 1, a, —a, 00} of the quadratic differentials with two double zeros of Section 13.21 are 
mapped to the poles in the z-plane {0, 1, ^1,^21 00}) and such that the critical graph (Figure H]) is 
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Figure 8: The projection of As on the ^i-plane. At various values of ^i, we show the section of the moduli space 
in the ^2-plane. The boxes all have the same scale. It is clearly visible that when we approach the boundaries Bi or 
B2, the section shrinks to zero. The thick boundaries of the sections are boundaries of the second kind, the others 
are of the first kind. 



mapped to a critical graph that is obtainable continuously from the graph of the right of Figure 
[2] by the vanishing of two lengths, and such that ^1 belongs to the region characterized by having 
positive imaginary part, real part smaller that 1/2 and absolute value greater than one. There are 
four maps obeying all these requirements. The first one is 



hi{w) 



a — w 
2a 



which maps 



hi{oo) = 00 , hi{—a) = 1 , hi{a) = 0, 
and ^1 and ^2 must therefore be the images of —1 and 1. 



ei = /ii(-i) = 

This corresponds to the curve Bi 



a + 1 
2a 



6 = /ii(l) 



a — 1 
2a 




where a{i) is given by (j3.37|) . The second map is 

a — Iw 



h2{w) 



1 



a + 1 w + 1 



(5.4) 
(5.5) 

(5.6) 
(5.7) 
(5.8) 
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It maps 



This defines the curve B2 



The third map is 



for which 



/i2(l) = 0, h2{-l) = oo, h2{-a) = l, 



6 = /i2(oo) 



a — 1 
a + l 



6 = /i2 



' ^ a{£) + 1 ' ^ 



"•2, 



2 — a 
a + 1 w — 1 ' 



(5.9) 

(5.10) 
(5.11) 



/i3(l) = oo, /i3(-l) = l, /13(a) = 
4a 

6 = /i3(-a) = , ,^9 > 6 = /i3(oo) 



(a + 1)2 



a + l 



We will call the corresponding curve C 



C 



4a{i) 
ia{£) + 1)2 ' 



" vr" 




0, - 


} 


L ' 2J 





(5.12) 



(5.13) 



Before we discuss the fourth map and the meaning of the third one, we need to look at what we 
have found so far. At the special point £ = 0, we recall (I3.38P that a(0) = 2 — a/S and thus we find 
that all three curves meet at the point r/3 (Figure [8]), where 



2^=0 



CI 



l + \/5 



1=0 



When £ = ^, the curve Bi meets the axis Re^i = ^ at the point 772- From (j3.38p we find 

1 32^/6 + 19^/15 

7)2 = 1 1 ■ 

' 2 54 

Also when £ = 5, the curve B2 meets the unit circle at the point 774, given by 



27+ (32\/6 + 19^/15) 
27- (32^/6 + 19^/15) 



(5.14) 



(5.15) 



(5.16) 



The situation for the curve C is a little bit different because for £ > £q k, 0.67622, the curve is 
outside of the allowed domain because it is inside the unit circle centered at the origin. But we 
can map the piece of curve with £ > £q back inside the domain by the map z ^ 1/z. So the curve 
C gives us actually two curves Ci and C2 (and the fourth map is just h4{w) = l/h^{w)). We have 
thus 



Ci 



Aa{£) 
(«W + 1) 



G [0,4] 



C2 



{a{£) + 1)^ 



4a(^) 



G 



vr 
''2 



(5.17) 
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At £ = io these two curves intersect on the unit circle at the point 

r/5 « -0.47019 + 0.88256 i. 
At last, when ^ = f , the curve C2 meets the axis Re^i = ^ at the point rji, with 

1 19^/15 



(5.18) 



(5.19) 



In Figure [9] we show the critical graphs of the quadratic differentials obtained by mapping the 





Figure 9: The critical graphs of the quadratic differentials obtained by mapping the quadratic differential of Figure 
|4]with, respectively, hi, /12 and /13. 



quadratic differential of Figure H] with, respectively, hi, /12 and h^. So they are respectively on Bi, 
B2, and €2- 

It remains now to understand the meaning of the curves Ci and C2. From Figure [9] we can be 
more specific by noting that on these curves, one quadrilateral and one triangle edge are zero. By 
exhausting all combinations on the prism of Figure [1] and Equ. (|1.5p . we see that this implies that 
we have also three triangle edges of length tt. We thus have an intersection of three boundaries of 
the first kind and one boundary of the second kind. This can happen if a boundary of the second 
kind of the section is reduced to a point. Therefore the curves Ci and C2 partition into 



three regions A 



(l,n) 



n 



1,2,3 in which the sections 5^^ have exactly n boundaries of the second 



kind. This partition is shown on Figure [HI where we can also verify the number of boundaries of 
the sections shown in the boxes. We do not prove that this partitioning is as shown, but we have 
checked it numerically beyond doubt. 



To end this section, we briefly sketch how we plot numerically the sections. We start by plotting 
the boundaries of the second kind; this is done with the quadratic differentials with one double zero, 
studied in Section 13.11 The end of the boundaries are detected when a length is vr, and accurately 
determined. We then start from one end (where, say, ii = tt) to plot a boundary of the first kind 
by finding the points ^2 for which ii = vr. We end either when another length (say ij) becomes vr, 
or when we end up very close to a boundary of the second kind. In the first case, we determine 
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accurately the corner (the intersection of two boundaries) and continue by looking at the points 
for which ij = tt. In the second case we go on by plotting the boundary of the first kind attached 
to the other end of the boundary of the second kind that we just met. We continue this process 
until we arrive at the other end of the second-kind boundary from which we started. The result 
will be a closed curve given by a set of points separated by an approximately fixed distance h; and 
the special points at the corners of the curve are always precisely given. 

6 The five-tachyon contact term 

To calculate the five-tachyon contact term, we take the general formula established in ^ for the 
A^-tachyon term (or use Equ. (|7.25p ). Namely 

2t/ _ (-1)^"' 2 f ^yf dx^Advi 1 

where = Xi + iyi. Since the five external states are the same, the integration over Vo,5 can be 
written as 120 times the integration over ^5. Therefore 

2 f - - 1 
K^Vt5 = ^ dxi Adyi Adx2 Ady2fi{(i,Ci,C2,C2) , 1^ = 1 , (6.2) 

It is now a good place to explain how we perform the numerical integration. There are five 
distinct steps in the process. 

Step 1 First we compute the boundary of the reduced moduli space. For this we draw a covering 
rectangular grid of (A^ + 1) x [N + 1) points on A^^\ the projection of on the ^i-plane. Then 
at each of the points ^1 of the grid we attempt to plot, in the ^2-plaiie, the boundary of the section 
S^^. In the ^2-plane, this boundary is represented as a list of points, and the space between two 
successive points is taken to be approximately h. It may happen that the algorithm fails to find a 
Strebel differential because the Newton method fails to converge (for example if we are close to a 
singular quadratic differential, near the boundaries Bi or ^62). In that case, the algorithm may be 
able to fix the problem by trying other initial values; if this still fails the consequence may be that 
the algorithm is unable to draw the section, but this is not dramatic, it simply keeps record of this 
failure and proceeds to the next point of the grid. We are using three different grids which have 
respectively {N,h) = (30,0.1), (70,0.05) and (100,0.03). On Figure lO] we show all the points of 
these three grids for which the section could be plotted. We see that we have some gaps near the 
boundary Bi, especially near the point 773 where it meets the boundary B2- This is not surprising 
because at this point the quadratic differential has two triple zeros, a very singular point indeed. 
But it is clear that as we increase the number of points on the grid, the gaps are becoming smaller. 
This is due to the fact that, going from one point to the next, the quadratic differential solver based 
on the Newton method, can use the solutions of the nearby points as seeds, and of course if they are 
closer the Newton method will have more chance to converge. The grid with N = 100 is however 
already at the limit of the computational power of a desktop computer. We notice also a few gaps 
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N = 100 



Figure 10: The projection of As covered with grids of (iV + 1) x (iV + 1) points with N = 30, 70 and 100. The dots 
represent the points of the grids for which the section could be computed. 

away from the boundaries, and most noticeably there seems to be a line of gaps in the middle of the 
grid A'' = 100. The failure there is actually not due to the Newton method but to the difficulty of 
plotting the boundary of a section when one very small boundary of the first kind is stuck between 
two other boundaries of the first kind. This could be fixed by improving the algorithm, but we 
don't really need to do this at this point because we will eventually use quadratic interpolation to 
fix these gaps. One may worry about the size of the gaps near Bi and especially near = 173, but 
we recall that near the boundaries Bi and B2 the sections become small, and therefore the gaps 
actually represent a small four-dimensional volume, and will be fixed by extrapolation. The error 
made by the interpolations and extrapolation will be estimated later. 

Step 2 Now we cover every section with a rectangular grid, with spacing s along the real and 
imaginary directions, and at every of these points we attempt to find the Strebel quadratic differen- 
tial. Again it is possible that some differentials are not found, but this will be fixed by interpolation 
(extrapolation). At every of these points, we store all the data of the quadratic differential, i.e. ^j, 
ttj, dtti/d^j, dai/d^j, i,j = 1, 2 and p/, / = 1, . . . , 5. For our three particular grids N = 30, 70 and 
100, the spacing s is chosen as respectively s = 0.1, 0.05 and 0.04. We also consider the same grid 
N = 70 with the smaller spacing s = 0.04. 

Step 3 We integrate on every sections. 



because the shapes of the sections are complicated, we use a Monte-Carlo technique and integrate 
the function 





(6.3) 




(6.4) 
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over a rectangle containing the section 5"^^. To evaluate at random points ^2 we use quadratic 
interpolation over the closest 3x3 sub-array of points of the grid. Actually, we find it more natural 
to interpolate log fi, and exponentiate the result of the interpolation. 



Step 4 We can finally perform the whole integration 

/ dxi Adyi Adx2 Ady2Ai(6,Ci,6,6) = / c^a^i Ad?/iM(^i,^"i). (6.5) 

We use again a Monte-Carlo technique, integrating 

M'(6,6) = (f^;'^"^).^^ ^^^^^'^ , (6.6) 
I U otherwise 

and integrate M' over a rectangle containing ^5^^ To evaluate M, we use again quadratic inter- 
polation of logM, and exponentiate the result of the interpolation; we observe that it is a better 
interpolation that if we had simply interpolated M itself. 

Step 5: Error estimation We would like to calculate an estimate of the numerical uncertainty 
on the result. There are four sources of error, namely 

1. The error due to the finite spacing h. In order to decide if ^2 £ •S'^n we check if ,^2 is inside 
the polygon formed by all points describing the boundary of 5*^^. This is basically a linear 
interpolation, so the error made will be of the order of /i^, a potentially large error. To reduce 
it, we are computing the integral twice, once in the normal way, and a second time with all 
the sections replaced with cruder ones with spacing 2h. Practically we are removing half 
the points of every sections but keep the corners (the intersections between boundaries). In 
order to isolate this source of errors from the other sources, it is important to use the same 
Monte-Carlo samples in both integrations (i.e. the sequence of random points must be the 
same). Let us denote K^V^s^h) and ^^145(2/1) the results of these two integrations, and let 
K^VfS be the exact value. Since the error is quadratic in h, we have 

kVj5(2/i) - kVj5 = 4 {K^Vt5(h) - kVjs) + 0{h^) , (6.7) 

or in other words 

kV^s = - {4:K^Vt5{h) - kVj5(2/i)) + 0{h^) , (6.8) 
3 

and it will turn out that we can neglect the residual error as it will be smaller than the other 
sources. 

2. The error done by the interpolation and extrapolation of fi. Let us write fj, = //inter + ^^i- We 
will then treat this source of error as systematic, in other words 



(72=/ dxi Adyi Adx2 /\dy2 \Sf^{£,i,^i,(,2,C2)\ ■ (6.9) 



A: 



5 



The quadratic interpolation is done with the routine polin2 of [19]. Based on Neville's 
algorithm, it can give an estimate of 6^{Ci,^i,S,2,^2)- 
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3. the error due to the interpolation and extrapolation of M. It is a systematic error as well. 
We write M = Minter + 5m and 

0-3 = dxiA dyi \5m\ (6.10) 

4. The error coming from the Monte-Carlo integrations. Those are purely statistical errors, and 
the errors coming from the integrations in Step 3 will cancel down to a negligible quantity 
in the final result if they are done with enough samples. We are thus only considering the 
Monte-Carlo error apic on the last integral of Step 4. 

So finally 

^'^Vs = 4 {4 + + <^Mc) , (6-11) 

TV 

and we will take enough samples in the Monte-Carlo integration so that the error will be dominated 
by the errors on interpolation and extrapolation, which depend on the finesse of the grid. 

We note that Steps 1 and 2 need to be done only once for every grid. Once these have been 
done, every integration requires only Steps 3, 4 and 5. 

We now state our results. We use four different grids, with (A^, h, s) respectively equal to 
(30,0.1,0.1), (70,0.05,0.05), (70,0.05,0.04) and (100,0.03,0.04). These last two grids are at the 
computational limit of a desktop computer with our C-|— |- code. Indeed it takes several days 
to compute aj(^i, ^1, ^2, '^2) and their derivatives and p/(Ci, ^1, ^2, '^2) at every points of the four- 
dimensional grid (70,0.05,0.04) (which has approximately 1.6 million points inside ^5, and fills 1.3 
GBytes of RAM). The grid (100,0.03,0.04) is actually too large to store the derivatives of Oj, so 
we can use it only for the computation of k^V^b. The computation of the five-dilaton contact term 
in Section [7] will have to be done with the grid (70,0.05,0.04). All the Monte-Carlo integrations 
are done with one iteration of the routine vegas of [19]. The integrations of Step 3 are done with 
10^ samples, and the integrations of Step 4 are done with 10^ samples. Our results are shown in 
Table [TJ It is reassuring that all four results are compatible within their error bounds. They are 



{N,h,s) 


(30,0.1,0.1) 


(70,0.05,0.05) 


(70,0.05,0.04) 


(100,0.03,0.04) 




10.0774 ± 0.065 


9.95537 ±0.013 


9.95197 ±0.010 


9.93412 ± 0.008 




10.0045 ± 0.065 


9.93535 ±0.013 


9.93252 ±0.010 


9.92648 ± 0.008 




9.980 ± 0.065 


9.929 ±0.013 


9.926 ± 0.010 


9.924 ±0.008 



Table 1: The results of the integrations and their uncertainties on our three different grids. 



not statistically independent however, and we should therefore not combine them. So we take the 
result from the finest grid as final answer 



9.924 ± 0.008 



(6.12) 
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7 The five-dilaton effective term 



The goal of this section is to check the vahdity of our computation. The quartic vertex computed 
in [7| was checked in [10] with marginal fields. A similar analysis is unfortunately not possible here 
because we can couple only an even number of marginal fields like a_iQ!_iCiCi|0) and the vanishing 
of the quintic term of its effective potential is trivial. However we can compute the amplitude of 
five dilatons \D), where 

ID) = (cic_i - cic„i) |0) , (7.1) 
and compare it to the prediction of the dilaton theorem. 

The dilaton theorem states that the effective potential of the dilaton should be identically zero. 
That the V^^ term vanishes is due to the fact that the cubic amplitude of three dilatons is zero 
because the ghost numbers do not work out. The vanishing of V^4 is nontrivial; it has been 
checked in [11] that the contribution from the quartic amplitude k^V^a cancels the contribution 
from Feynman diagrams with two cubic vertices. This cancellation was found to be good to about 
0.2%, furnishing a very good evidence that the quartic computations are done right. In this section, 
we want to check the vanishing of V^^ . There are two contributions to this effective term shown 
on Figure fTTj namely the contact term /t^V^s, and the Feynman term C^s 




Figure 11: The term d!^ of the dilaton effective potential. The fields {^i) that propagate in the second diagram, 
are the tachyon and all the massive fields, but not the dilaton. 

kV/ = kV,5+Q5. (7.2) 

We will start by evaluating the Feynman contribution, then we will compute the contact term with 
the machinery developed in this paper. 

7.1 The Feynman contribution 

The three-string vertex can couple only an even number of states with asymmetric ghost numbers 
like the dilaton. So the diagram consisting of three three-string vertices doesn't contribute, and 
the only Feynman diagram contributing to C^^s is therefore the one shown on the right of Figure 
nn To get the right relative sign between this diagram and the contact term, we must keep track 
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of all the coefficients. Since the vertices bring in the amplitude a factor —i, and the propagator a 
factor i, we have 

H)5!Crf5 = (f\ Y,i-i){D,D,D,^,}{{} {-M-')^^{-i){D,D,^j} , 

^^'i' = -Y2 ^ {D,D,D,^,}{M-')^^{D,D,^,} , (7.3) 

where (— M^^) is the zero-momentum propagator for all the fields except the dilaton. M is therefore 
given by the quadratic term 

Mi, = {^i\cQQB\^j) , i,j^2. (7.4) 

Note that we work in a non-orthonormal basis, and M is therefore not diagonal. Qs must be 
evaluated with level truncation, we will do explicitly the first few levels. First it is convenient to 
decompose into the sum of contributions at every level 

Q3(£) = ^c,5(0. (7.5) 

£'<£ 

And we will write Mi for the matrix M with only fields of level i, and we define the row vectors 

Ce = {{D,D,^f^}) 

Q, ^ ({AAA^O), (7.6) 

where the index i runs over all fields of level i. With these notations, we can write 

cA^ = -\^Qi-M-^.Cj. (7.7) 

At level zero, only the tachyon propagates, and we thus have 

Q5(0) = A ATjMo-HA AT} = -l(-5.716)(-^)(-^) = 0.4020. (7.8) 

At level four, the four fields are (we are using the notations of [15j ) 

1^3) = (6-2ClC_2Cl - 6_2CiC__2Cl) |0) , |*4) = C_iC_i|0) 

1*5) = L_2CiL_2Ci|0) , l^-e) = (c_iZ_2Ci - c_iL_2Ci) |0) . (7.9) 



And we have 



M4 = diag(-4,2,338, -52) 



48' 432 ' 72^ 

^4 = (-4.547,-1.811,-22.33,10.69) (7.10) 
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The first two lines are done with standard techniques, and the third hne has been calculated in 
[15] . we refer the reader to [15] and [11] for the details of this computation. With these values ()7.7p 
gives the contribution of level four 

c^5 (4) = -0.03996 (7.11) 



which gives the total contribution 

Q5(4) =Q5(0) + Crf5 (4) = 0.3620. 

At level six, the fields are 



(7.12) 



^7) = (6_2C_2C_iCi - 6_2C-2C_lCi) |0) , j^s) = (-^^-2C_3Cl - L_2C„3Ci) |0) 
^'g) = L^2L-2 (C-lCi - C_iCi) |0) , l^'io) = (6_3CiC_3Ci - l-sCiC-sCi) |0) 

= (6-3ClL_2C-lCi - 6_3ClL_2C-lCl) |0) , 1^12) = C_2C-2|0) , |^'13) = L_3ClZ_3Cl 1 0) 
^14) = 6_2C-lCi6_2C_lCi|0) , 1^15) = (c_2-^^-3Cl - C-2L-3C1) |0) 

l^'ie) = (c_26-2C-lCi - C_26-2C-lCi) |0) , I^'it) = (L_3Ci6_2C_iCi - L_3Ci6_2C_iCi) |0) , 

and we have 



Me 

And we find 



(8) 



/ -104 \ 

-1352 

0-8 

V -104 / 



4 
10816 
4 



416 
8 
416 



100 
729 



(-8,0,0,0,0,1,0,4,0,4,0) 



(7.13) 



(-2.103, 2.363, -20.92, -0.4412, 7.090, 0.1251, -49.67, 0.4998, 4.325, 1.603, -8.649) 



c^5 (6) = -0.03606 , 

which gives the total contribution 

Q5(6) =Q5(4) + Crf5(6) = 0.3259. 



(7.14) 



(7.15) 



With the techniques developed in [3 [HI [TTl [TOl US] , we can evaluate {D,D,D,^i} for fields 
of level up to ten. The results for C^s {£), £ = 0, 4, 6, 8, 10 are shown in Table [2j We will follow the 



i 





4 


6 


8 


10 


CM 


0.4020 


0.3620 


0.3259 


0.3167 


0.3112 



Table 2: The value of with propagating fields up to level 



method used in [T^ to extrapolate the above results to ^ = oo. Namely we use a fit of the form 

/: 



cAi) = k + j^- 



(7.16) 
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In fitting C^i{i!.), the authors of [TT] found that 7 = 3 give the best results. Here by fitting the four 
values Q5(£), I = 4,6,8, 10 with (|7.16p . we find that 7 = 2.51. We take this as an evidence that 
the data should be fitted with 

/i 



Q3W = /o + ^. (7.17) 



With this last fit we find 



Qs =^5(00) = 0.3060 . (7.18) 



7.2 The contact term 

The multilinear A^-string function at genus zero is given by [U [9] 

{^'i,...,^,v}= (^^J dAiA...AdA2(iv-3)(S|K^Aj...K^A2(^_3))l*i)---l^^)- (7-19) 

The Xi are 2(A^ — 3) real coordinates of the reduced moduli space Vq^n of the A^-punctured spheres, 
and the surface state (S| corresponds to the A^-punctured sphere given by the parameters Aj. The 
antighost insertions h{v\.) are given by 



N 00 



b{v..) = Y. E {HU€'+HlJg)) , where Hm = f ^^J^^'^ ■ (7.20) 

7=lm=-l I * 

The hi(wi; Ai, . . . A2(Ar_3)) are the maps from the local coordinates wi to the sphere, which are 
going to be described a little later. Now we follow [11] and rewrite the antighost insertions in a 
more convenient form. First we rename the coordinates of the moduli space as 

ii = Xi + iyi , where Xi = X2i-i , Vi = X2i , « = 1, . . . , A^ - 3 , (7.21) 

where the are the complex coordinates naturally used in the quadratic differentials formalism. 
We have 

dxi A dyi b{v^jb{vyj = d^A diBiB* = -2i dxi A % BiB* , (7.22) 

where 

N 00 N 00 

7=1 m=—l 1=1 m=—l 

and the coefficients Bf and C/ are given in terms of derivatives with respect to and 



jji _ f dw 1 1 dhi (-,1 _ I dw 1 1 dhi 
And we arrive at the formula 

{^-1, ...,^n} = I dxiAdyiA...A dx^^s A dyiv-3(S| (BB*)^ . . . {BB*)j^_^ \^n) ■ 

(7.25) 
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We will now consider the case of interest, N = 5 and \^i) = \D), where \D) is the dilaton (j7.ip . 
we will thus only need the coefficients and C-^ with m = —1 and m = 1. To evaluate them, 
we need to expand the maps hj to order 



z = hi{wi;^i,(i,^2, 6) = zi + piwi + p](ii wj + phi + 0{wj) , (7.26) 

where all the coefficients on the right-hand side depend on ^i, ^i, ^2 and ^2- The zj are the positions 
of the finite poles on the z-plane, namely 

zi = , ^2 = f , zz = ii , 2:4 = 6 • (7.27) 

For the puncture at infinity, we must use the coordinate t = 1/z 

t = h5{w5; 6, e'l, ^2, 6) =P5W5 + pI/35 wI + pi75 wl + OiwI) . (7.28) 

We can now use (|7.24p to express the coefficients that we need. We find 

Bi-i = — h,i , B2 I = — 6ii , C- I = 
oil 2 ' d42 2 



e/ = - 67/ . (7.30) 



where 



The coefficients /3/ and 7/ can be calculated from the expression of the quadratic differential ([2 
and (I2.12p . First we expand Lp = (j){z){dz)'^ around the puncture zi 

and knowing that in the local coordinates (p takes the canonical form 



1 



we find the relations 



^ = (7.32) 

Pi = \^-i ' li = YQi^{rUf + ^{ri?)- (7.33) 



And after expanding (t){z) to obtain and Tq, we find 

/33 
/34 



266 

ai + «2 + 26 - + 26 - ei 

2(6 - i)(6 - 1) 
Qi-26'-^i + 6(3 + a2 + 26) 

26(6-i)(6-6) 

ai - 6^ + 266 + 6(3 + a2 - 26) 



26(6-i)(6-6) 

/35 = -^(2 + a2 + 6 + 6) , (7.34) 
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and for ej = 8/3] — 67/, we find 

_ -5{u - ai)2 + 12 ((^3 + ^2^(1 + -ai{s + t)+ ^^(l - + ^1) + *(-3 - 02 - 6 + CD) 



€2 



8^2 

-5{u - 2s - ai - 02)^ - 12 (5s -7u + 2w + t{4s - 5 - u - t) - 02(5 - 2) - ai(2s - t - 3)) 

8(t - s + 1)2 

^3 = iTZTT A77 TTW (4e?(106 - 15 - a2 - 146) - ai(5ai + 16^? + 2^(3 + 5a2 - 26) 

ol4ua - 6jia - J-jj 

+2(6 - 56)6) + ei(12(l + 6) - Sei) - 2t(12 + 96 - 5a26 + 2^1) 
+e?(15 - 18a2 - 5oi + 246 - 8026 + 44^1)) 

^4 = WTT, rVc TTW + 10^? - 126 + 4t + 2(-5a2 - 3 - 86)6) 

8(6(6 - 6j(6 - ijj 

+26'(-6 + (9 - 5a2)6 - 22^1) + 8t(3 - (3 - 02)6 + 7^1) 
-C|(15 - 18a2 - 5ai - 606 - 4a26 + + fii^^i - 12 + 46)) 

^5 = ^ (16 + 12ai + is + 7u + 2t + a2(-5a2 - 8 + 2s)) , (7.35) 
8 

where we have defined 

s = 6 + 6, t = 66, ^ = 6' + e2', w = ei+f2- (7.36) 

We will also need to compute a few correlators. To see which ones, it is best to start expanding 
the antighost insertions. It is easily seen from the expression for the dilaton and the fact that only 
the correlators with ghost number (3, 3) do not vanish, that wc need either three 6's and one b, or 
one b and three fe's; and as already mentioned we only need antighost oscillators with indices — 1 
or 1. Let us expand {BB*)^ {BB*)2 in the following way: wc take all terms in the first factor that 
have either two 6's or two 5's and mixed terms in the second factor. The rest of the expression, 
with mixed terms in the first factor, can be found simply by changing some indices 

{BB*),{BB% = 

1=1 1=1 Ijtj I^J 



Bi-ib^l E Biib? - ^,b^l E Kib? + E MiA%'^ + (1 - 2 , ^ - ^) , 

7^4 7^4 IjtJ J 

(7.37) 



where 



Ml' ^Bl,Bi,-Cl^CJ, (7.38) 



and the last term is obtained by changing as indicated the left subscripts of B and C, the subscript 
of M, and all superscripts. Noting that 

b.i\D) = c_i|0) , bi\D) = -ci|0) , 5_i|D) = -c_i|0) , bi\D) = ci|0) , (7.39) 
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we see from (j7.37p . that we need the following open correlators 

Aij ^ {{c^iCi)^'\c^^l)o , Bij ^ ((c_iCi)(^),4^))o (7.40) 
CuK = {cf^c^\cf\ , D,jK ^ {c%c['\cf\ , E,jK ^ {c%c%c[^\. 

The conventions for the closed correlators are as in 

{c{zi)c{z2)c{zz)c{wi)c{w2)c{wz)) = -2(c(zi )c(z2 )c(2;3))o {c{wi)c{w2)c{wz)) o , (7.41) 
and the open correlator is 

{c{zi)c{z2)c{z^))o = {Zi - Z2) {Zi - Z3) {Z2 - z^) . (7.42) 

Ajj and Bjj were already calculated in [llj. For the other correlators, we can use either the 
conformal transformation of the operator ^5^c(z) corresponding to c„i (as in [H]), or conservation 
laws as in [15j. We find when I , J, K ^ 5 (with the definition Zjj = zi — zj): 

- 2f3if3jzij + ^ejzij{l - Pizi, 



(7.43) 



Au = 


PJ - 


Bij = 


1 

PJ 


ClJK = 


1 


z 

PIPJPK 


DijK = 




PJPK V 


ElJK = 


'''' u 

PK V 



- PliziK + Zlj)zjK + ^ejzjjz, 

Pj + PiPj{zij + ziK - zjk) + ^ejzjK - ^Piej{zij + zik)zjk+ 
+^e/ {^ziK + Pjzik{zjk - zij) + ^ejzjjziKZjK^ ^ , I,J,K ^b. 

The cases when oscillators are at infinity must be treated separately. Since i?f _i and Cf_^ are 
zero, we cannot have a c_i alone at infinity. We therefore need 

= pj(]^ej{P^ + zj)- Pj 



Bl5 
ClJ5 



^ , B,j = -{zj + P,) 

P5 PJ 
ZJI 

PIPJP5 



DlJ5 = (Pl - \^IZIJ 

PJP5 V 2 

Eij5 = ^ (\Piej - \eiPj - \eiejzu^ , (7.44) 
P5 \2 2 4 J 

and we note that Cjjk is totally antisymmetric, Djjk is antisymmetric in J and K, and Ejjk is 
antisymmetric in I and J. 
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We can now write the integrand of (j7.25p with five dilatons from (j7.37p and the definitions (j7.40p . 
{n{BB%{BB%\D)\D)\D)\D)\D) = 



4Re < 



{Ak3 {D^ijCi 



+Ak4 (-Ds/jC/ ii?2,i ~ BjjCf iB2^i + BjiCliBl ]) 
+Wj [-Eui [cliBl, - C^iBl^ + Bl, {AnCl, - Ai^Cl,) + Bl, {AmCI, - AnCl,))] 

+- Yl {^(-^3/jC{iM/^ + BijCl.M^'' - BjjCl.Ml'' 

{~^3IjBIiC2^i + -B/ji?|iC^i - -Bj/i?2,iC|i)^ - CijkAl3BIiC2^iB^i^ 

+- V \bI]^ {-DiijCi iM/^ + BijCl iM/^ - BjiCi^Mf'' 

Pa ^ — ' L ^ ' ' ' 

+^2,1 {~^4:ijBl iC^i + BjjBf^^C^i - BjiB{^^Cf i)^ - CijkAl4B{iC(iB2^^ 



+ CijkBtl {BliCl.M^'^ + Bl^Cl^M^"^) 



(7.45) 



With (17.29p . (j7.34p . (|7.35p . (|7.43p and (|7.44p . this expression can be expressed in terms of ^j, Oj, 
dtti/d^j, dai/d£,j and p/, and we can therefore integrate it numericahy over the reduced moduh 
space of five-punctured spheres. Since we have five times the same state, the result is simply 120 
times the integral on ^5, thus 

kV^s = ^ {D, D, D,D,D} = ^ [ dxi A dyi A dx2 A dy2(S| {BB*\ {BB*)2 \D)\D)\D)\D)\D) . 

(7.46) 

We do the integration as in Section [6l on three different grids. The results are shown in Table [3l 
They are compatible within their error bounds, and we take again the result from the finest grid 



[N,h,s) 


(30,0.1,0.1) 


(70,0.05,0.05) 


(70,0.05,0.04) 


K^Vd5{2h) 


-0.30666 ± 0.0093 


-0.30768 ± 0.0023 


-0.30759 ± 0.0016 




-0.30343 ± 0.0093 


-0.30666 ± 0.0023 


-0.30660 ± 0.0016 




-0.3024 ± 0.0093 


-0.3063 ± 0.0023 


-0.3063 ±0.0016 



Table 3: The results of the integration (|7.46p and its uncertainty on three different grids. 



as final answer 



-0.3063 ±0.0016 



(7.47) 



We see that ()7.18p and (j7.47p cancel each other with a precision of about 0.1%, well within the 
error bound of 0.5% on the contact term (|7.47p . This is solid evidence that our computations are 
reliable. 
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8 Conclusions and prospects 



In the light of the verification successfully made in Section [TJ that the effective potential of the 
dilaton vanishes at order five, we can claim that the techniques described in this paper work well, 
that the reduced moduli space is understood and described right, and that our implementation of 
the algorithm gives reliable results. In particular we trust the value (|6.12p obtained for the contact 
term of five tachyons. It is also a good check of the consistency of CSFT itself. 

We are able to estimate the uncertainty made in the computation of terms (see (|6.12p and (|7.47p ). 
However we want to mention here that we have been quite conservative in this estimation. In 
particular the errors on the fits were treated as systematic instead of independent. This is only half 
right because we expect interpolations made from different samples to be more or less independent, 
although two interpolations at two nearby points using the same sample are not independent. The 
overestimation of the error seems confirmed by the values in Tables [J and [3l However one should 
be careful in this respect because, as we already mentioned, we don't think that the errors made 
on different grids are independent either. So unless we find a better way of estimating the error, we 
will stick to our conservative estimation. Any way, we emphasize that the errors on any given term 
will probably be no more than 0.5%, the error we find for the five-dilaton term (which is calculated 
from a very long expression). And this precision is probably enough to do level truncation. 

The effect of the five-tachyon term ()6.12p on the stable vacuum is studied in [15]. However this 
term of level zero isn't enough to draw any conclusion. For this, it will be necessary to compute 
other terms at higher level. At level two there is only one term, namely four tachyons and one 
dilaton ¥^4^. At level four we'll have five terms, Vfi^2 and V^i^. where i = 3, . . . ,6, are the 
four scalar fields at level four. The main difficulty in computing these terms will be that when 
we have different external states, the integrations on the 120 pieces of reduced moduli space won't 
be all equal, and it will require some (straightforward but lengthy) work to express them as one 
integral over ^5; the extreme case being when we have five different states, we will have 60 different 
integrals (complex conjugation is always trivial and divides the number of integrations by two). 
But the computation of interactions to level four is certainly not that bad. By automatizing the 
computation of oscillator algebra and correlators, as in [15], it might be possible to compute higher 
levels as well. 

At last, it would be useful to make our numerical data available, so that readers can use it to 
make their own computations of quintic terms. Ideally we would like to create fits of the boundaries 
of the moduli space and the parameters of the quadratic differentials ai and 02 and the mapping 
radii, that could hold in a paper and be entered in a computer in a reasonable amount of time. This 
is under study, but we are a bit pessimistic given the dimensionality and rather complicated shape 
of the reduced moduli space (our data describing its boundary is about 600 megabytes large). 
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